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We study quantum and classical Hanbury Brown- Twiss correlations in waveguide lattices. We 
develop a theory for the propagation of photon pairs in the lattice, predicting the emergence of 
nontrivial quantum interferences unique to lattice systems. Experimentally, we observe the clas- 
sical counterpart of these interferences using intensity correlation measurements. We discuss the 
correspondence between the classical and quantum correlations, and consider path-entangled input 
states which do not have a classical analogue. Our results demonstrate that waveguide lattices can 
be used as a robust and highly controllable tool for manipulating quantum states, and offer new 
ways of studying the quantum properties of light. 
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Lattices of coupled waveguides have been shown to be 
extremely versatile in manipulating the flow of light [l|, [2] . 
Such systems enabled direct observation of optical ana- 
logues of many fundamental quantum mechanical effects 
such as Bloch oscillations [J, |4j], Anderson localization 
[5|, |y, L0|, quantum Zeno effect [8|], quantum walks [9| 
and many others [l|, [2|, llO|- However, these effects were 
all related to the wave properties of light rather than 
to its particle nature. To observe the quantum prop- 
erties of light, one should consider correlations between 
single photons. Here we show that photon pairs propa- 
gating in waveguide lattices develop nontrivial quantum 
correlations unique to lattice systems. We experimen- 
tally observe some of these features, albeit with reduced 
contrast, using Hanbury Brown- Twiss intensity correla- 
tion measurements with phase- averaged coherent states. 
Our results suggest that waveguide lattices can be used 
as a robust and highly controllable tool for manipulating 
quantum states in an integrated manner, offering new 
ways of studying quantum properties of light in struc- 
tured media. 

Consider a lattice that is made of identical optical 
waveguides, each supporting a single transverse mode, 
as shown schematically in Fig. la. The evolution of the 
quantized electromagnetic field in each waveguide is given 
by the Heisenberg equation for the bosonic creation op- 
erator a) . For a lattice with nearest-neighbors couplings, 
the Heisenberg equation for the creation operator of the 
k th waveguide is given by: 
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where z is the spatial coordinate along the propaga- 
tions axis, (3 is the propagation constant of the waveg- 
uides, Ck,k±i are the coupling constants between adja- 
cent waveguides, and c/n is the speed of light in the 
medium. The creation operators at any point along the 
propagation are calculated by integrating Eq. (pQ): 
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Uk,i(z) is a unitary transformation given by calculat- 
ing the exponent of the coupling matrix izCk,i, which de- 
scribes the amplitude for the transition of a single pho- 
ton from waveguide I to waveguide k. Since any input 
state can be expressed with the creation operators a] 
and the vacuum state |0), the evolution of non classical 
states along the lattice can be calculated using Eq. ([2]). 
When a single photon is coupled to waveguide Z, the in- 
put state a] |0) = \l) l will evolve to the superposition 
^Z k Uj" k |l) fe , where \l) k is the state of a single photon 
occupying waveguide k. However, measurements of the 
probability distribution of single photons are not enough 
to reveal the quantum properties of light, as the prob- 
ability distribution of a single photon \Uk,i(z)\ 2 evolves 
in the same way as the intensity distribution of classical 
light [Si, ll2|, [l3| . The quantum mechanical properties of 
light are observed when correlations between the propa- 
gating photons are considered. In this Letter we focus on 
the evolution of the photon-number correlation function 
r g?r = (a^alardq) , when two indistinguishable photons 
are injected into the lattice. 

We start by analyzing the simplest example of two 
coupled waveguides. The coupling matrix in this case 

and the transformation Uk,i(z) is 
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isin(Cz) cos(Cz) 
pier acts as a beam splitter, with the reflection and 
transmission coefficients varying continuously along the 
propagation. If two photons are injected to the cou- 
pler, one to each waveguide, the average photon number 
at each of the waveguides is constant since 711(2) (z) = 

\ a i(2) a i(2)) = l^ii 1 2 + W\2\ 2 = 1. The nonclassical na- 
ture of the light is revealed by considering Ti^, the prob- 
ability to detect exactly one photon at each waveguide 
(a coincidence measurement). Using Eq. ([2]) we obtain 
ri, 2 (*) = \U11U22 + U 12 U 2 i\ 2 = Cos 2 (2Cz). Since two 
paths lead to the final state of one photon at each waveg- 
uide, they interfere and the probability for a coincidence 
measurement oscillates along the propagation. After 




Figure 1: (a) A schematic view of the waveguide lattice used 
in the experiment. The red arrow represents the input light 
beam, (b) The calculated probability distribution (n q (z)) of 
a single photon injected to the central waveguide of a periodic 
lattice, as a function of the propagation distance. The photon 
couples coherently from each waveguide to its neighbors, and 
the probability distribution concentrates at two outer lobes. 
(c) The calculated correlation matrices r g , r , representing the 
probability to detect at the output of the lattice exactly one 
photon at waveguide r and one photon at waveguide g, when 
both photons are coupled to a single waveguide at the center 
of the lattice, i.e. |<^o) = l/v^o 2 |0). This matrix is a simple 
product of two single-photon distributions, thus showing no 
quantum interference. The grey bars are obtained by sum- 
ming the matrix along one axis, representing the results of a 
single photon measurement. 



propagating exactly half a coupling length zbs = k/4C, 

we find Ti^(zbs) = 0. At this point the two photons 

bunch, and are found together in either one of the two 

< — > 
waveguides. U (zbs) is identical to the transformation 

of a symmetric beam splitter, and the coincidence mea- 
surement vanishes in the same manner as in the Hong- 
Ou-Mandel (HOM) interferometer [14]. Since a pair of 
coupled waveguides is equivalent to a beam splitter, it is 
possible to cascade several of them in order to implement 
quantum gates in an integrated manner [15|. A lattice 
of many coupled waveguides enriches the variety of cor- 
relations obtained in integrated structures, as we show 
bellow. 

We now turn to study the quantum properties of a pe- 
riodic lattice with a large number of identical waveguides, 
where all the coupling constants are equal C n>n ±i = C. 
As long as the photons are far from the boundaries of the 



lattice, Eq. ([2j) yields U q ,k(z) 



*J q -k(2Cz), where J q 



is the q th Bessel function [2|, ll6|]. When a single photon 
is coupled to waveguide fc, it will evolve to waveguide q 
with a probability n q = \U q ^(z)\ 2 = J q _ k (2Cz) 2 . The 
photon spreads across the lattice by coupling from one 
waveguide to its neighbors in a pattern characterized by 
two strong 'ballistic' lobes, as shown in Fig. lb. If a sec- 
ond photon is coupled to another waveguide Z, then the 
average photon number at waveguide q is simply the in- 
coherent sum n q =< a q a q >= J q -k(2Cz) 2 + J q -i(2Cz) 2 . 
Once again, the quantum nature of light is revealed by 
considering the correlations between the two photons. In 
the following we study the correlation T Qjr for three dis- 
tinct two-photon input states: i) Both photons are cou- 
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Figure 2: Quantum and classical correlations in waveguide 
lattices, (a) The correlation matrix r g , r when the photons 
are coupled to two adjacent waveguides i.e. \ipi) = %a\ |0). 
The two photons exhibit bunching, and will emerge from the 
same side of the lattice, (b) The correlation matrix when 
the two photons are coupled to two waveguides separated by 
one waveguide, ^2) = al 1 a{ |0) . Here the two photons will 
emerge either both from the lobes, or both from the center. 
(c) Measured classical intensity correlations T^r, correspond- 
ing to the \(fi) input state, (d) Measured classical intensity 
correlations f^f, corresponding to the ^2) input state. 



pled to a single waveguide at the center of the lattice, i.e. 
\<po) = l/y/2a,Q |0) (ii) The two photons are coupled to 
two adjacent waveguides |y?i) = aja J |0) and in) The two 
photons are coupled to two waveguides, separated by one 
waveguide, ^2) = a_ r al |0). The correlation matrix T qjr 
represents the probability to detect one photon at waveg- 
uide q and its twin photon at waveguide r/g. The prob- 
ability to detect both photons at the same waveguide q 
is given by T qiq /2. 

Fig. lc depicts the correlation matrix Y q ^ r at the out- 
put of the lattice, when the two photons are coupled to 
the same input waveguide (the |(/?o) input state). In this 
case, there is no interference and the correlation matrix 
is just a product of the two classical probability distribu- 
tions, T q ^ r = 2\U q oU r o\ 2 - The correlation map is charac- 
terized by four strong lobes at the corners of the matrix, 
resulting from the tendency of the photons to propagate 
in the ballistic directions. 

When the two photons are coupled to two neighbor- 
ing sites, i.e. the \tpi) state, the correlation map changes 
considerably as shown in Fig. 2a. The most obvious fea- 
ture is the vanishing of the two 'off- diagonal' lobes: the 
photons tend to bunch to the same lobe. This can be 
thought of as a generalized HOM interference. Two paths 
lead to a coincidence measurement between waveguide q 
and waveguide r: either the photon from waveguide 
propagates to waveguide q and the photon from waveg- 
uide 1 to waveguide r, or vice versa - from waveguide 



to waveguide r and from waveguide 1 to waveguide 
q [17|]. These paths are complex and involve hopping 
of the photons between many waveguides, nevertheless 
they interfere and the correlation matrix is thus given by 
r g , r = \U q oU r i + U q iU r o\ 2 - The destructive interference 
which leads to vanishing of the off-diagonal lobes can be 
mathematically traced to the inherent 90° phase shift as- 
sociated with nearest-neighbors coupling. We note that 
since the photons tunnel between the waveguides con- 
tinuously, the visibility of the quantum interference only 
weakly depends on the overall length of the lattice. 

The four lobes are recovered when the photons are ini- 
tiated in the ^2) state, i.e. with one waveguide sepa- 
ration between the input sites (Fig 2b). However, this 
state also contains strong non-classical features - note 
the differences between Fig. lc and 2b. The photon pair 
exhibit bunching but with a different symmetry: if one 
photon is detected in between the lobes, the probabil- 
ity to detect the second photon in a lobe vanishes, even 
though a single photon is most likely to reach the lobes. 
Similarly, if one photon is detected in a lobe, it is certain 
that the other photon is also in a lobe. It is important to 
note that in the above examples, the quantum interfer- 
ence emerges since the two photons are indistinguishable. 
Interferences with distinguishable particles are possible, 
but require an entangled input state [181 ]. 

Some of the special features of these quantum me- 
chanical correlations can be captured using intensity cor- 
relation measurements with classical light [13, 120, l21 |. 
The quantum mechanical probability to detect one pho- 
ton at waveguide q and its twin photon at waveguide r 
is related to the classical intensity correlation function 
Tq/ r = (I q I r ) where (•) denotes statistical (or tempo- 
ral) averaging. This non-local intensity correlation func- 
tion is usually discussed in the context of the classical 
Hanbury Brown- Twiss effect and its quantum interpre- 
tation [ll|, [22| . We studied experimentally intensity cor- 
relations at the output of a periodic waveguide lattice. 
The lattice of 89 identical waveguides was fabricated on 
an AlGaAs substrate using e-beam lithography, followed 
by reactive ion etching. Each waveguide is 8mm long 
and 4/xra wide [23j]. The tunneling parameter between 
sites C is determined by the etch depth of the sample 
(1.3/xra) and by the distance between neighboring waveg- 
uides (4/im) and was measured to be 290ra -1 . We used 
an OPO (Spectra-Physics, OPAL) pumped by a mode- 
locked Ti:Sapphire laser (Spectra-Physics Tsunami) to 
generate 150/s pulses, at a wavelength of 1530nm with 
80MHz repetition rate. The average power was on the 
order of 0.1m W, thus nonlinear effects were negligible. 

Two-photon input states can be mimicked by injecting 
into the lattice two phase- averaged coherent states: two 
coherent states with the same mean photon number 
and a fluctuating relative phase [13, [25] . Thus, we 
injected two coherent beams into two different sites of 
the lattice [24|], and randomized their relative phase with 




Figure 3: Quantum correlation maps r g , r for path-entangled 
input states. (a) The two photons are injected to- 



gether to either of two neighboring waveguides, lip 



(+)) 



\ (°i +°o ) |0)- The correlation is significant only in the 
'off-diagonal' peaks, which indicates that the two photons 
will emerge from opposite sides of the lattice, (b) The cor- 
relation map for an input state with two photons in either 
of two waveguides with one waveguide separation, where 

there is a 7r-phase between the two possibilities, w^ 



±(a?-al\)\0). 



a spatial light modulator. For each phase realization, 
the intensity profile at the output facet of the slab 
was imaged on an infrared camera, and the intensity 
correlations between all the waveguides were computed. 

(c) 

The intensity-correlation function T q ^ was obtained by 
averaging the measured correlations over many phase 
realizations. The measured intensity correlations for 
the nearest neighbors (mimicking the \ip±) state) and 
next-nearest neighbors (^2) state) inputs are presented 
in figures 2c and 2d, correspondingly. The patterns are 
strikingly similar to the corresponding quantum corre- 
lations T qjr , except for the reduced contrast. Indeed, 
for two incoherent sources coupled to waveguides k and 
I 7^ fc, the classical correlation is given by T q c : l(z) = 
I§ {\UqiU rk + U qk Uri\ 2 + \U q iU r i\ 2 + | U qk U rk | 2 ) , where 
io is the intensity coupled to each waveguide. The 
last two terms are responsible for the reduced contrast. 
Thus, while the 'off-diagonal' lobes practically vanish for 
the quantum input state \ipi), a straightforward calcu- 
lation shows that for classical light T q c Jr > ^ y ^q,q^r,r , 
in agreement with our experimental results. Still 
these classical Hanbury Brown- Twiss type intensity 
correlations do echo many of the special features of the 
quantum correlations of Fig. 2a, b. 

The great potential of waveguide lattices for quantum 
information probably lies in the extension of these con- 
cepts to non-uniform lattices in order to specially de- 
sign the correlation properties, and in the utilization 
of these multiport systems for handling more complex 
quantum states. As an example for the latter, consider 
the propagation of a path-entangled input state with 
two photons in either of two neighboring waveguides, 
|V^ + )) = \ la{ + aj ) |0). The calculated correlation 
map for this case is presented in Fig 3a. The corre- 



lation in the 'diagonal' peaks completely vanishes and 
is significant only in the 'off- diagonal' peaks. Accord- 
ingly, the two photons will always separate and emerge 
from different sides of the lattice. The corresponding 
correlation map violates the Cauchy-Schwarz inequality 
Tq^r < y / Tg^TV~r , indicating that this case has no appar- 
ent classical analog. As another illustration, consider the 

state l?/^ - ^) = \ (a[ — a_A |0). Now the photons are 
in either of two next-nearest neighboring waveguides, yet 
with a 7r-phase shift. The correlation map shown in Fig. 
3b reveals that in this case, one photon will always reach 
a lobe while the other will always reach the center. 

In this Letter we studied the evolution of photon pairs 
along periodic lattices, and have shown that the result- 
ing quantum correlations strongly depend on the input 
states. We compared our predictions with a classical 
wave theory, and experimentally demonstrated that some 
features can be obtained using classical intensity correla- 
tion measurements, yet with reduced contrast. The cor- 
respondence between the classical and quantum nature 
of light can be further studied by considering the evolu- 
tion of quantum correlations in the presence of dephas- 
ing, which can be introduced via lattice inhomogeneities. 
Furthermore, waveguide lattices offer new possibilities as 
they can be designed in ways that are not feasible using 
bulk or fiber optical systems. It will be especially inter- 
esting to study the effect of such lattices on the propaga- 
tion of other types of non-classical light, such as squeezed 
states [26] and cat states. 
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